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Abstract 

The Higgs potential of 2IIDM keeps its generic form under the group of transfor- 
mation GL(2,C), which is larger than the usually considered reparametrization group 
SU{2). This reparametrization symmetry induces the Minkowski space structure in the 
orbit space of 2IIDM. Exploiting this property, we present a geometric analysis of the 
number and properties of stationary points of the most general 2IIDM potential. In 
particular, we prove that charge-breaking and neutral vacua never coexist in 2HDM and 
establish conditions when the most general explicitly CP-conserving Higgs potential has 
spontaneously CP-violating minima. We also define the prototypical model of a given 
2HDM, which has six free parameters less than the original one but still contains all the 
essential physics. Our analysis avoids manipulation with high-order algebraic equations. 



1 Introduction 

The Standard Model relies on the Higgs mechanism of the electroweak symmetry breaking. 
Its simplest realization is based on a single weak isodoublet of scalar fields, which couple to 
the gauge and matter fields and self- interact via the quartic potential, for review see [2j . 
Extended versions of the Higgs mechanisms are based on more elaborate scalar sectors. The 
two-Higgs-doublet model [3], where one introduces two Higgs doublets 0i and 02, is one of 
the most economic extensions of the Higgs sector beyond the Standard Model. This model 
has been extensively studied in literature from various points of view, see [H HI El [6] and 
references therein. The minimal supersymmetric extension of the Standard Model (MSSM) 
uses precisely a specific version of the 2HDM to break the electroweak symmetry, 
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The Higgs potential of the most general 2HDM Vh = V2 + V4 is conventionally parametrized 

as 
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with 14 free parameters: real m1^, m22, Ai, A2, A3, A4 and complex m^2? A5, Ag, A7. The large 
number of free parameters makes the analysis of the most general 2HDM and its phenomeno- 
logical consequences rather complicated. It suffices to say that even the equations for the 
stationary points and the mass matrices of the most general 2HDM are too cumbersome to 
solve them explicitly or to extract any useful information in a straightforward way. Appendix 
A of [1] serves as a remarkable illustration of this point. 

Fortunately, in many phenomenological applications one does not actually need to consider 
the most general 2HDM. For example, the spontaneous CP- violation takes place in the Higgs 
sector of 2HDM even with = and real A5, Ag, Ay, as originally studied in [5]. In an even 
simpler case m^2 = A5 = Ag = A7 = there is still much room for very rich phenomenology 
since the two doublets can couple to fermions in different ways. The MSSM is also based upon 
a very specific version of 2HDM with Ai = A2 = — 2A3 and A5 = Ag = A7 = 0. In all these cases 
the straightforward algebra is usually sufficient for the complete treatment of the EWSB. 

Still, there are issues that go beyond these simple versions of 2HDM and call upon the 
study of the most general 2HDM. First, one might be interested in understanding what can 
and what cannot happen, in principle, if the Higgs mechanism involves two doublets. Next, 
when discussing natural sets of the free parameters of the Higgs potential, one runs into the 
question of crucial and, in some sense, redundant parameters. In other words, one wants to 
know the "space of different physical situations" realized in 2HDM, not the space of different 
free parameters. Attempting to link different-looking Higgs potentials with the same physical 
content with some reparametrization transformation (see below), one unavoidably plunges into 
the analysis of the most general 2HDM. Finally, one can go even further and consider more 
evolved Higgs sectors, like 3HDM, where the straightforward algebra strongly limits one's 
insight into phenomenological possibilities. 

The problem that emerges is not to calculate precisely the vacuum expectation values 
(v.e.v.'s) of the fields, the Higgs masses etc., but to understand the general structure in the 
space of all 2HDM's. 



1.1 Reparametrization invariance 

For a long time, the knowledge of this structure had been very limited. One would set 
some parameters to zero, and occasionally observe that some of the remaining ones were 
"more important" than the others. The situation began to change several years ago with the 
observation of the reparametrization freedom and the development of the basis-independent 
methods, [HI E]- 

In general, a theory has reparametrization invariance if two different sets of parameters 
lead to the same observables. In 2HDM one has the Higgs potential ([1]) that depends on Higgs 
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fields (pa and coupling constants Aj, mfj, V = V{\i, mf-] (pa). The physical observables are 
functions of these coupling constants, but not fields themselves. The key observation is that 
if one performs a reparametrization of the potential, 

Aj, mj Xi, , (2) 
then the resulting Higgs potential V{\i, mf^-; (pa), via a transformation of the Higgs fields 

(pa ^ (pa = Ra/3(pf3 , 

can in some cases be brought back to the form with original coefficients V{Xi, mfj] (pa)- When 
this happens, the essential physics encoded in V^(Ai, mf^; (pa) and V{\i, mf^; (pa) is the same. 
In other words, if one performs transformation ([2]) and, simultaneously, inverse transformation 
of the fields, (pa — > {R~^)a/3(pi3, one obtains exactly the same expression for the potential as 
before. 

By construction, the Higgs fields and the coupling constants live in different spaces. It 
is convenient to think of the Higgs potential as a scalar product between the "tensors" com- 
posed of fields and "tensors" composed of coupling constants. The initial form of potential 
corresponds to the explicit coordinate representation of these tensors in some specific basis. 
Rotating the basis, one induces the above transformations of the fields and of the coordinates, 
which leaves the potential and the physics invariant. 

The reparametrization symmetry offers the freedom to choose the basis for description of 
the Higgs potential. The quantities that are invariant under basis rotations reflect the physics, 
while quantities that depend on the basis transformations contain redundant information 
about the approach chosen. 

The standard way to implement the reparametrization freedom is to consider the dou- 
blets (pi and (p2 as up and down components of a "hyperspinor" $ and perform a unitary 
transformation $ —>■ f/$, U G U{2). All the 2HDM's whose potentials are related by uni- 
tary transformations from U{2) correspond to the same physical theory written in different 
bases. Obviously, for the basis-independent analysis, one must use the most general form of 
the potential. 

Representing the Higgs potential as 

V = + iz,,,,($t$,)($t$,) , (3) 

with tensors Yab and Zatcd being the short notation for the mass terms and quartic couplings, 
respectively, one can construct from them many new tensors. The complete traces of these 
tensors are the SU{2) invariants. The physical observables must be constructed from these 
invariants. The parameters that do change under the SU{2) transformation are redundant; 
they depend on the way we describe the model but not on the model itself. A heavy machinery 
can be developed along these lines, see [6] for details. 

1.2 Discovering structures in the space of two-Higgs-doublet mod- 
els 

Representation ([3]), although compact, sheds little light on the structure of the space of all 
2HDM and, therefore, on the properties of a generic model. For example, a great calculational 
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effort was undertaken in [8j to establisli tlie complete set of basis-independent conditions 
for tfie explicit CP-conservation in the most general 2HDM. This analysis revealed four real 
invariants, but the meaning and origin of them remained a mystery. 

The structure of the space of all 2HDM became clearer after [10], where a simple yet 
efficient group-theoretic analysis of the 2HDM was performed. Tensors Yab and Zabcd were 
decomposed into irreducible representations of SU{2) and then were mapped onto the irre- 
ducible representations of the group 5*0(3). The complicated tensorial analysis was reduced 
to the simple linear algebra in 3D linear vector space with one real symmetric matrix and 
two real vectors. Basis invariant conditions for the hidden Z2 symmetry and for the explicit 
CP-conservation could be easily derived in these linear- algebraic terms. 

Recently, a similar approach was used in [Tl] to extend such geometric analysis to the 
models with Higgs doublets. From the very beginning, the author of |Tl] switches from the 
fundamental to the adjoint representation of SU{2), or SU{N), and derives the criteria for the 
absence of explicit and spontaneous CP- violation in these models. Another realization of such 
analysis was presented of |12j . Authors of this work went on to find and study the stationary 
points of the 2HDM potential using interesting algebraic techniques. The authors of another 
recent paper [T3] also studied a generic NHDM, aiming at discovering general properties and 
the underlying structure of this model. 

It is interesting to note that back in 1970-1980's there was much activity on mathematical 
properties of various realizations of the symmetry breaking Higgs mechanisms. It was under- 
stood that the problem of minimization of some group- invariant potential is simplified if one 
switches from the space of Higgs fields to the orbit space, [H]. This idea was exploited in [15j 
and in |T6j to study the minima of a Higgs potential invariant under the Lie group G with 
Higgs fields transforming under various representations of this group. In electroweak theory, 
C = SU(2) X U{1)ew^ but this simple case has never been analyzed in detail under the most 
general circumstances. For a recent application of these old ideas to 2HDM see [T7] . 

The recent analyses as well as the older mathematical work make it clear that in spite of 
the messy parametrization ([1]), the Higgs potential in 2HDM must have a very simple internal 
structure. In this paper we make a further step to uncover and put to work this structure. 

In Section 2 we note that the Higgs potential ([T]) keeps its generic form under the action 
of a larger group of transformations between the two doublets, GL(2, C). This leads us to the 
group 50(1, 3) of the transformation of the adjoint representation, and induces the Minkowski 
space structure in the orbit space of 2HDM. Even more important, it induces the Minkowski 
space structure just in the space of 2HDM's themselves. We conclude Section 2 by defining 
a prototypical model of a given generic 2HDM, which involves six parameters less but still 
contains all the essential physics. 

We then use this formalism to study the existence and the number of the stationary 
points that correspond to the charge-breaking (Section 3) and neutral (Section 4) vacua, as 
well as the mass matrices at these points. We do this in geometric terms, avoiding explicit 
manipulation with high-order algebraic equations. Then, in Section 5 we derive a criterion 
when the most general 2HDM with explicit CP-conservation has spontaneously CP-violating 
solutions. Finally, in Section 6 we draw conclusions. Three Appendices provide additional 
details. 
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2 Minkowski space of 2HDM 



2.1 Extending the reparametrization group 

The Higgs potential ([1]), defined in the 8- dimensional space of Higgs fields, depends on them 
through their electroweak-invariant combinations (EW-scalars) Xij = {(f)l(j)j), i,j = 1,2. The 
points in the space of Higgs fields that have the same values of Xij can be related by a unitary- 
transformation from the electroweak group SU{2) x U{1)ew and therefore form an orbit. 
Points within each orbit are indistinguishable for the Higgs potential, and one can think of 
the Higgs potential as defined in the orbit space. 

Since the two Higgs doublets 0i and 02 in 2HDM have the same quantum numbers, we treat 
them as the up and down components of a single contravariant hyperspinor $ = 02)"^, 
which realizes the fundamental representation of SU{2). The quartic part of the potential 
can be viewed as a twice-covariant, twice-contravariant tensor, with coefficients Aj being its 
coordinates in a chosen basis. The rotation from one basis to another is realized by a unitary 
transformation from SU (2), which induces the corresponding transformation both of the Higgs 
fields $ ^ $' (and, consequently, the basis tensors) and of the coordinates, A A', leaving the 
physics invariant, see details in ^U\. This reparametrization group offers freedom in describing 
the Higgs potential and, in general, the entire 2HDM. 

It is very convenient to switch from the fundamental to adjoint representation of SU{2), 
As 2 ® 2 = 3 © 1, the quantities Xij form a singlet and a triplet. Using the well-known 
SU{2) 5*0(3) mapping, one maps them into a scalar and a real-valued vector: 

/ i4<P^) + (0l02) \ 
ro = ($t$) = (0l0i) + , r, = (<|.ta,$) = -^[(^I^a) - (0^0i)] • (4) 

V (0l0i) - (0^02) / 

The Higgs potential can be now written as 

V = -MiVi - Moro + AjTiVj + BiTiVo + Crl , (5) 

where Mj and Mq contain the mass coefficients, while Ajj, Bi and C are composed of quartic 
couplings Aj, see details in [TUl [HI [12] • 

Representation ([5]) already displays some structure in the space of all possible 2HDM, i.e. 
the space of all free parameters of the potential. An SU{2) rotation induces a corresponding 
5*0(3) rotation of the basis in 3D space, under which r,, Mj, Bi transform as vectors, Aij 
transforms as a symmetric tensor, but the value of the potential remains the same. The phys- 
ically relevant quantities, the scalars of 50(3), are not just these vectors and tensors but their 
scalar products and traces. 



We now make a step forward and consider the largest group of invertible linear transfor- 
mations GL{2, C) of a complex- valued 2- vector. Since, for instance, the quartic part of the 
Higgs potential contains all possible fourth order terms, an arbitrary linear transformation 
between the two doublets keeps it intact, only up to reparametrization. Therefore, the group 
under which the Higgs potential is reparametrization-invariant is GL{2,C), not just U{2) as 
usually assumed. 

The subgroup C* (i.e. overall multiplication by a non-zero complex number) of GL(2,C), 
generated by a unit operator, is redundant for the description of the potential. Multiphcation 
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of all the fields by the same real non-zero constant leads to rescaling of all the observables, 
without changing the structure of the model, while the global phase rotations have no effect 
on the potential. It is its factorgroup SL{2, C) that embraces all non-trivial transformations 
and generates interesting symmetries. 

For the adjoint representation, SL{2, C) induces the proper Lorentz group 5*0(1, 3). Apart 
from the 3D rotations, induced by SU{2), one also has boosts along the three axes. The scalar 
Tq and vector rj now become parts of a single irreducible representation of S0{1, 3): 



(6) 



Thus, the orbit space, which is given by all possible four-vectors r^, is equipped with the 
Minkowski space structure. The covariant and contravariant vectors are related, as usual, 
by g^lv = diag(l, — 1, — 1, — 1). For convenience, when describing geometry of this space, we 
will coin the standard space-time terms "future" , "past" , "time-like" and "space-like" , whose 
definitions are obvious. 

Now, the Higgs potential in the orbit space can be written in a very compact form: 



(7) 



where 
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A1+A2 1 \ 
— V A3 

-Re (Ae + A7) 

Im (Ae + A7) 



-Re (As + A7) Im (Ag + A7) 

A4 + Re A5 — Im A5 

— Im A5 A4 — Re A5 

Re (Ae - A7) -Im (Ag - A7) 



Ai — A2 



Re (Ae - A7) 
-Im (Ae - A7) 
A3 



A1+A2 
2 



This way of presenting the potential appeared in [TT]; however index /i was assumed to be just 
a short notation for a pair "singlet + triplet" . The wider group of spinor transformation was 
not exploited. 



2.2 Extended reparametrization group and kinetic term 

The extended reparametrization group SL{2, C) modifies the Higgs kinetic term. Introducing 

p/^ = (9„<I>)V'^(9"$) , (9) 

where a indicates the true space-time coordinates, one can rewrite the kinetic term as 

i^ = V. (10) 

In the standard form, = (1, 0, 0, 0), which sets a preferred "reference frame" in the space 
of p'^ and, therefore, in the orbit space. Upon an SL{2, C) transformation of the Higgs fields, 
spacelike components of Kfj_ become non-zero, and the canonical form of the kinetic part is 
broken. 

The structure of the kinetic term does not affect the number of the stationary points of 
the potential, however it does affect the position of these stationary points and the matrices of 
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second derivatives calculated around them. If ip^ are the neutral scalar fields, then the second 
derivative matrix 

1 9V 

M„b = - , a = l,...,8, (11) 

represents the true mass matrix of the Higgs bosons only when calculated in the above preferred 
"reference frame". However, the true mass matrix and the second derivative matrix in any 
Lorentz frame share a common property: the positive definiteness. 

Proposition 1. Any SL{2,C) transformation of the Higgs fields preserves the positive (semi)- 
definiteness of the second derivative matrix. 

Proof. Consider each of the complex Higgs fields as a pair of real Higgs fields. An SL{2, C) 
transformation induces an <S'L(4, R) transformation of the resulting eight real Higgs fields 

The mass matrix is then transformed as 

Suppose M'^fj is positive (semi) definite; then M^^f^qaqb > (or > 0) for any vector g^- Denote 
Qa = qa'(R~^)a'a- Siuce map Ra'a IS surjectivc, then QaQbMab > (or > 0) for any Qa- 
Therefore, Mab is also positive (semi)definite. □ 

This proposition states nothing but the fact that a local minimum remains a local minimum 
under any rotation and stretching of the coordinates in the space of Higgs fields. 

Thus, the fact that kinetic term is not conserved under SL(2, C) transformations has little 
relevance to the study of the overall properties of 2HDM. The formalism to be developed here 
is particularly useful to study the following issues: 

1. the number of stationary points of the potential and their charge-breaking/neutral na- 
ture; 

2. the minimum/saddle point classification of the stationary points and the number of the 
Goldstone modes; 

3. spontaneous CP-violation in an exphcitly CP-conserving 2HDM; 

4. possible phase transitions under continuous variation of the parameters of the potential. 

In the present paper we pay attention to the first three issues and only occasionally comment 
on possible patterns of phase transitions. 

2.3 The orbit space 

The orbit space in 2HDM is not the entire Minkowski space. The square of the 4-vector 
is invariant under any proper Lorentz transformation and, due to the Schwartz lemma, is 
non-negative 



2 u 2 2 

r = r'^r^ = Tq - r ^ 



{<t>\<t>l){<j>\<t>2) - (0l02)(0^0l)l > . (12) 
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Therefore, the physical field configurations lie inside and on the border of the future light- 
cone (LC^) in the Minkowski space. The cone shape of the orbit space in 2HDM was known 
before, see for example [17] , however, its meaning in the context of Minkowski space was never 
realized. 

The surface of = 0, corresponds to the situation when the two Higgs doublets 

01 and 02 are proportional to each other. In particular, if a vector indicates a stationary 
point (the vacuum) of the potential, then = means that the vacuum is electrically neutral. 

> can be realized only when the two doublets are not proportional; a vacuum solution 
with > would correspond to the charge-breaking vacuum. 



2.4 Properties of A^"" 

The positivity constraint on the Higgs potential requires it to be bounded from below. These 
constraints are usually presented as a list of inequalities among different A's, see, for example, 
[5]. In our language the positivity constraint is given by a single statement: 

the tensor h.^^^ is positive definite on the future light-cone , (13) 

which, as is shown in Appendix |Al is equivalent to the following set of requirements: 

• tensor A.^" is diagonalizable by an 5*0(1,3) transformation, 

• the timelike eigenvalue Aq is positive, 

• all spacelike eigenvalues Aj are smaller than Aq. 

In other words, there always exists an 5*0(1,3) transformation that brings A^'' to 



with Aq > and Aq > Ai, A2, A3 . 



Note that Aj are bounded only from above. The negative values of Aj with arbitrary large 
absolute values are allowed. In particular, all A, < if and only if A.^'^ is positive definite in 
the entire Minkowski space. 

Below we focus only on the physically relevant A^'^, which satisfy (fT3l) . 
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2.5 The prototypical model 



Since is diagonalizable, one can always perform such a 50(1, 3) rotation that brings it to 
the diagonal form. By doing this, one arrives at the prototypical model of the initial 2HDM 
model. This model is defined by the diagonal K^'^ , or explicitly, by 



Ai — A2 — A , A, 

The potential of the prototypical model 
1 



A7 = , Im A5 
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has only four free parameters in the quartic potential A, A3, A4, A5 and a generic set of the 
mass terms. 

Thanks to Proposition [1], issues 1-4 listed in Section 12.21 depend on Lorentz-invariant 
quantities, and, therefore, can be studied at the level of prototypical model. The answers 
to these questions are the sam^ for the entire six-dimensional manifold of 2HDM's that are 
linked to the given prototypical model by S0{1, 3) transformation^. Therefore, if one intends 
to study issues 1-4 in the most general 2HDM, one can immediately get rid off Xq, A7, and 
Im A5 and set Ai = A2, without any lack of generalit}{§. 

We stress that we do not write down explicitly the expressions relating Aj with the coupling 
constants Aj of the generic 2HDM. This would require solving a fourth-order equation, and the 
expressions would be just impractical. The point is that in order to understand the coarse- 
grained structure of a generic 2HDM we do not need this exact transformation law. We do 
not calculate nor use the numerical values of vacuum expectation values or masses of physical 
Higgs bosons. If one is interested in properties of a 2HDM with particular values of coupling 
constants, one can always diagonalize A'^'^ numerically and immediately check what class of 
models the given model belongs to. 



3 Charge-breaking vacuum 
3.1 Stationary points 

Representation ([7]) is a convenient starting point to study the properties of the vacuum in 
2HDM. The orbit space one obtains upon the spin map 0j — >■ r'^ is a region with border {LC~^) 
in the Minkowski space. Two possibilities must be considered in search for stationary points 
(or better to say, stationary orbits) of the potential. The first one is when a stationary orbit 
in the Higgs field space corresponds also to a stationary point in the orbit space. In a generic 
case it lies inside LC^ and corresponds to the charge-breaking vacuum. The second possibility 
is that a stationary orbit in the space of Higgs fields corresponds to a non- stationary point on 
the border of the orbit space, i.e. on LC"*", which corresponds to a neutral vacuum. In this 
Section we consider the former case, while Section 4 will be devoted to the latter one. 

The stationary orbit is composed of stationary points, and the conditions for the stationary 
point of the Higgs potential are 

^ = ^^ = <C. = 0, .= 1,2, where < = , e = -M^ + k^^r^. (16) 

Eq. ( fT6l) is a system of four equations since both real and imaginary parts are set to zero. 
Introducing the light-cone vectors n'^ = (1,0,0,1), = (1,0,0,-1) and the "transverse" 
unit vectors = (0, 1, 0, 0), 63 = (0, 0, 1, 0), one can rewrite (fT6|) as 

0i<CM + 02(ei'-ze^)CM = O, 02ri'lCM + </'i(e5' + ^e^)C = O. (17) 



^ There is a subtlety with spontaneous CP- violation; see Section [5] 

^One can further reduce the number of essential free parameters without losing the generality by rescaling 
the fields and switching from A, A3, A4, A5 to 1, A3/A, A4/A, A5/A. 

^In principle, the boost transformations do not remove the free parameters but rather move them to the 
kinetic sector of the Higgs lagrangian. However, as was shown before, the exact structure of the kinetic term 
has little relevance to the issues we study here. This comment is due to Roman Li. 
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Consider the case of charge breaking vacuum. The two doublets are non-zero and are not 
proportional to each other, so that each of f|T71) sphts into a pair of conditions that correspond 
to vanishing coefficients in front of (pi and 02 separately. One obtains 

CX = 0, C«^ = 0, C,e>^ = 0, Ce^ = 0, (18) 

from which one gets = 0, or 

A^''r^ = M^. (19) 

This is an inhomogeneous system of linear equations. In order to classify all situations here, 
let us introduce the image of LC^ under the action of tensor K^^^: 

LC' = {x^\x^ = A^^r^ for all r'^ G LC+} . (20) 

The following criterion holds for existence of the charge- breaking stationary point (s) in the 
orbit space: 

Proposition 2. Charge-breaking stationary points exist if and only if lies inside the 
lightcone image LC . 

Proof. The proof follows from Eq. ( fT9l) and from the properties of the operator A.^'^ described 
in Appendix El Indeed, a point lies inside a cone if and only if the image of this point lies 
inside the image of the cone. Since lies inside LC~^, K^^Ty lies inside LC", and in order for 
physical solutions of f|T9l) to exist, must lie inside LC . □ 

The number of solutions depends on whether K^^^ is invertible. 

1. If all eigenvalues of K.^" are non-zero, then is an invertible operator. The solution 
of (fl9l) exists, is unique and is given by 

= (A-l)A'^M^ (21) 

It corresponds to a physical solution only when lies inside LC . The value of the 
potential at this point is 

V = -]^{X~yM^My. 

In the language of the prototypical model (i.e. in the frame where K^'^ is diagonal), one 
has 

Since lies inside LC , this value is always negative. 

2. If at least one of Aj is zero, then K.^^ is not invertible. The kernel of A^*^, Ker A, is a 
subspace in the orbit space, and LC collapses to a manifold of smaller dimension. If 

^ LC", there are no solutions. Otherwise, one has a continuum of charged solutions 
whose dimension is equal to the dimension of Ker A. If one consider and M^, 
restrictions of and to the space orthogonal to Ker A, then K^" is invertible and 
a generic solution can be represented as 

r^' = {K-^ Y^'My + p^' , where G Ker A . (22) 

The value of the potential at these points is 
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3.2 Mass matrix 



The stationary points found above can correspond either to the minima or the saddle points 
of the potential The non-trivial maxima cannot exist for a polynomial potential of fourth 
degree bounded from below. The simplest proof is given by the directional minima technique 
discussed in [15] . One first chooses some direction in the Higgs field space and considers a ray 
starting from the origin along this direction. The Higgs potential depends on the coordinate r 
along this ray as = Mr^ + Ar^ with some real M and some strictly positive A. The nontrivial 
stationary point of this function can be only a minimum; this is the directional minimum. 

Now one considers all possible directions in the Higgs field space and the manifold of 
all such directional minima. The true stationary point of the potential must belong to this 
manifold; therefore at this stationary point there always exists at least one direction along 
which the potential is concave upwards. Therefore a non-trivial maximum of the potential, 
which requires that the potential be strictly concave down at this point, is impossible. 

In order to find the nature of each stationary point, one needs to calculate the second 
derivative matrix. Upon differentiating in respect to the fields, the electroweak index "opens 
up", and one must trace it down carefully. The Higgs fields will be now written as 
with i = 1,2 labeling the first and second doublet, while a = 1,2 labeling the up and down 
components of each doublet. It is also convenient to switch explicitly to real and imaginary 
parts of the fields. Since = for the charged vacuum stationary points, one gets: 

-- 4:1m df^Jm dip Af,^, 

- ARed'^Jmdlf.A^,. (23) 

Now, to merge these equations into one, consider the real vector (pi^a with i = 1, . . . , 4 defined 
as: 

(fi^a = (Re Im Re 02,a, Im (p2,a) ■ (24) 

The second derivative matrices fl25D are then 





dRe 


-pi^a <9Re 






d^V 




dim 


i>i,a dim 






d^V 




dRe 


i>i,a dim 





where 



(9V 










ei 




\ 







ri\_ 










ei 




n_ 
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^2 







M 
?2_ 


/ 



(26) 



As a final step, we merge all y^j^Q, into a single 8D real vector 

V^a = (v^i,!' V^i.s) , a = 1, . . . , 8 , (27) 
and for any pair of indices a = {i, a} and a' = {i', a'} define 

Ka, = Ef, ■ . (28) 
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Then the second derivative matrix takes form 



(29) 



Note that S^^, is invariant under rotations that mix the real and imaginary parts of the 
same field or mix the electroweak indices a and (5. This leads to the observation that matrix 
(129|) is covariant under electroweak transformations, which is, of course, a plain consequence 
of the scalar potential V{(j)) being an EW-scalar. This fact means that the second derivative 
matrix must have four zero modes, which correspond to simultaneous EW-rotations inside 
doublets 01 and 02- 

In order to study the positive-definiteness of the matrix fl25|) . we exploit Proposition [1] and 
perform a convenient 5*0(1,3) rotation. Since the charged vacuum solution r'^ lies inside the 



LC^, one can always perform such a 5*0(1, 3) rotation that makes 



0, 0, 0), with real 



u. Then, one is free to choose a representing point on this orbit by performing a convenient 
EW rotation. A particularly simple representing point is 



V2 V « 



u 



The second derivative matrix fl29l) written in this reference frame is 



2u^ 





(ei|ei) 


-(ei|e2) 


(ei|n_ 


) 





(ei|^4 


-) 





(ei|ei) 


(ei|e2) 


\ 




-(e2|ei) 


(e2|e2) 


-{e2\n 


-) 





-{e2\n 


+) 





-(e2|ei) 


-(e2|e2) 






(n_|ei) 


-(n_|e2) 


{n^\n- 


-) 





(n_ ri_ 







(n_|ei) 


(?2_|e2) 






































(n+|ei) 






-) 












(n+|ei) 


(n+\e2) 






































(ei|ei) 


-(ei|e2) 


(ei|n_ 


) 












(ei|ei) 


(ei|e2) 




V 


(e2|ei) 


-(e2|e2) 


(62!^- 


) 







) 





(e2|ei) 


(e2|e2) 


/ 



where we introduced short notation 



(30) 



(31) 



Upon obvious rotations inside four pair of coordinates, one arrives at the above mentioned 
four Goldstone modes, which in this particular basis correspond to 



!)i,2 ± Im 02,i) 



^1,1 + Im 02,2) , 



1 



^1,1 



Re 



^2,2 j 



When gauge sector is taken into account, these modes will become the longitudinal components 
of the four gauge bosons, including the photon. The other four modes have the following second 
derivative matrix 



Mf. 



1 d^V 

2 difidipj 



( 



2u^ 



(eo|eo) 
(ei|eo) 
(e2|eo) 
(ealeo) 



(eo|ei) 
(ei|ei) 
(e2|ei) 
(eslei) 



(eo|e2) 
(ei|e2) 
(e2|e2) 
(e3|e2) 



(eoles) 
(eilea) 
(e2|e3) 
(ealea) 



(32) 
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with 

< = i« + n^_) = (1, 0, 0, 0) , = i« - n^_) = (0, 0, 0, 1) , 

and e^, 63 defined at the beginning of this section. We remind that the matrix (l32l) is not the 
mass matrix of the physical Higgs bosons, see discussion in Sect. 12.21 but is hnked to it by a 
certain 5*17(4, R) transformation. 

Despite the very suggestive form of fl32|) . there is no simple relation between the sets of 
eigenvalues of A^'^ and the mass matrix. Indeed, the matrix in fl32|) has form 

[ ^00 —^oi 

while the matrix A'V = Qu'v is 

/ Aoo -Ao. \ 
V -A., ) ■ 

However, it is easy to establish the following general condition for the charge-breaking mini- 
mum: 

Proposition 3. The charge-breaking stationary point is a minimum if and only if is 
positive definite in the entire Minkowski space. 

Proof. Suppose Pi = {po, pi, p2, Ps) is a normahzed eigenvector of (132|) with the eigenvalue 
M2. Then 

= Mf^piPj = u'^A^'P^P, , where P^ = poe^ + pic^ + p^e'i + pae^ . 

If A^^ is positive definite in the entire space, K^^P^Py > VP^, and all the eigenvalues of the 
mass matrix are positive. The stationary point is, therefore, a minimum, and according to 
Proposition [1], it remains a minimum when transformed back to the "preferred frame" with 
canonical kinetic term. 

Inversely, if A^'^ is not positive definite, then a vector exists such that A'^'^Q^Q,^ < 0. 
Since {cg, e^, 62, 63} form a basis in the Minkowski space, can be uniquely represented as 

= Qoeo + qie^ + ^26^ + ^363- Then, 

u'A'^'-Q^Q, = M^^q.qj. 

Representing qi in the basis of eigenvectors {pf"^} of the mass matrix, gj = C(a)Pj^°-*, one 
simplifies this to 

a 

which can be satisfied only when at least one of mass matrix eigenvalues is negative. □ 
This Proposition might be also guessed from the expression for the second derivative matrix 

Thus, in order for the charge-breaking stationary point to be a minimum, all Aj must be 
non-positive. In the language of the prototypical model, the charge-breaking minimum takes 
place, when 

A3<A, A4>0, |A5|<A4. (33) 

Note, finally, that if A^^" has at least one trivial eigenvalue, then the mass matrix will also 
have trivial eigenvalues, i.e. additional Goldstone modes. This can be most easily seen from 
the fact that detA^^ = -detM^-. 
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4 Neutral vacuum 



4.1 Strategy for finding solutions and general expression for the 
mass matrix 

The neutral vacuum solutions are located on the border of the orbit space. Therefore a neutral 
stationary orbit of the potential in the space of Higgs fields does not necessarily correspond 
to a stationary point in the orbit space. However, in order to keep the terminology simple, we 
will use term "neutral stationary point" for "point in the orbit space that corresponds to the 
neutral stationary orbit in the space of Higgs fields" . 

Since LC^ is invariant under 5*0(1, 3), one simplifies the analysis from the start by making 
an appropriate boost to align the "time-like" principal axis of K^'^ with the future line, so that 
all = 0. After this, one will need to perform only 3D rotations. 

For the neutral vacuum solution, the two Higgs doublets are proportional, and equations 
( fT6l) are not reduced anymore into four equations ( |T8|) . In order to find all the solutions, we 



use the following strategy. Take a point = |f^(l, n) lying on LC~^ and perform a 3D 
rotation that makes r'^ = ^v^n'^, which corresponds to setting 02 = 0. In this way one obtains 
three instead of four equations (ITSl) : the only difference is that now Cm'^- ^ ^ = C^+- 
The solution of the resulting equation 

- A^^ ■ v^n+ ^ - C< = . (34) 

exists, and is unique, if lies on the half plane A^'^ ■ — (n'^, parameterized by real 

numbers > and (. Note that this plane never collapses to a line, because A'^" does not 
have any non-trivial eigenvectors on the LC~^. The values of f ^ and ( in this Lorentz frame 
are 

where we used short notation ( !3T1) . Since + nt = 2eQ and due to A''* = 0, one can simplify 
C to 

C = -^MX + ^^'Ao - ^v^{n+\n+) = ^v^A^ - Mq , (36) 

so that if one fixes Mq, the variables ( and are linearly dependent. The depth of the 
potential at the stationary point is 



V = —yM.nl = ■ (37) 



1 2,. M 

-V MuU'l = ; — ^-—^ , 

4 2 (n+|n+) 

In contrast to the values of and (, this depth is a Lorentz-independent quantity. 

Repeat now the same check for every possible lightcone vector. Geometrically, we take 
every point = ^v'^{l, n), find its image (r')^ under the action of A^'^, attach a line generated 
by at this point, and check if lies on it. As we continuously check all possible lightcone 
vectors, the half-plane generated by f ^ and ( sweeps some region in the orbit space. Therefore, 
the problem of finding all stationary point of the Higgs potential is reformulated as the problem 
of finding this region and establishing how many times each point is swept as we check all 
possible vectors G LC^. 

This problem will be solved below in geometric terms, see Proposition [7l But one fact is 
plain to see now: 
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Proposition 4. Non-trivial solutions to the stationary point problem can exist only if lies 
outside the past light-cone LC~ . 

Proof. The representation of v"^ in (!35l) requires that M^n^ > for some lightcone vector n'^. 
But if M'^ lies in LC~ , then for any vector G LC~^ the scalar product M^n'^ is negative 
(this becomes obvious after making a boost to align with the past direction). Therefore, 
if lies inside the past light-cone LC~ , there are no nontrivial solutions. □ 

The second derivative matrix at a generic neutral stationary point can be also found. In 
the case of the neutral vacuum solution, 7^ 0. Therefore, in the second derivative of the 
potential one needs to differentiate not only C^, but also df. One obtains 





dRe 


i>i^a dR^ 


hp 




d^V 




dim 


■pi,a dim 






9V 




dRe 


■pi,a dim 





+ 4Redf^^Redlf,A^,, 



2Re (al^-C 
2Re (cr'^jCf,) 6^13 + 4Im d'^ Jm d^j^^ Af.^ , 
-2Im «^.C^) + 4Re <Jm d'<,A 



(38) 



Again, one can merge these matrices into a single expression using fields (pa, a 
defined in ([27D and E^^^, defined in ([28D: 

d'V 



and 



(^n^, and then 



difa difb 

To simplify the analysis, perform a 3D rotatioii^ to make 
an EW rotation to choose the representative point 

of the orbit. This choice leads to the following expression for the second derivative matrix 



(39) 



d^V 
difa dipb 



/o 




















\ 
































4C 


























4C 































2f d) 


2f2(ra+ 62) 






































2f ^(d n+) 





4C + 2t;2(ei|ei) 


4( + 2v^{e2\e2) J 


Vo 











2v^(e2|n+) 





2v^{e2\ei) 



(40) 



One sees that the modes corresponding to charged {a = P = 1) and neutral {a = (3 = 2) 
excitations split. For the charged modes, one has two Goldstones and two massive modes with 
eigenvalues 4^. Among the four neutral modes there is one Goldstone, Im 0i 2, while the other 
three modes, ipi = (Re 01^2, Re 02,2, 02,2), have the following mass matrix 



1 d^V 
= -- — — 

2 dcpiOcpj 




{n+\n+) (n+ld) (?2+|e2) 
+ I (ei|n+) (ei|ei) (ei|e2) 
(e2|n+) (e2|ei) (e2|e2) 



(41) 



*Note that A'^ , in general, is not diagonalized by this rotation. 



15 



Since (n_|_|n+) is positive, there is always a positive eigenvalue of M^. With = and 
the mass matrix simplifies to 

/Ao+(e3|e3) (ealei) (eglea) \ /O \ 

Mij = v^l (eilea) Ao + (ei|ei) (dlea) + -2Afo . (42) 

V (e2|e3) (eslei) Ao+(e2|e2)/ VO -2Mo / 

In order for the stationary point to be minimum, one must simultaneously have C > ^ind 
all eigenvalues of (H2!) positive. Although the eigenvalues of both matrices in (H2l) are known, 
v^(Ao — Aj) > for the first matrix and 0, —2Mq for the second matrix, the eigenvalues of 
their sum are not related in a simple way to them, since these matrices do not commute. The 
explicit expressions for the eigenvalues can be written down, but they do not give any simple 
minimum criterion in the general case. However, when A^'^ is positive definite in the entire 
Minkowski space, the analysis simplifies and the following important fact can be easily proved: 

Proposition 5. The charge-breaking and neutral minima never coexist in 2HDM. 

Proof. According to Proposition [2], a charge-breaking solution exists if and only if lies 
inside lightcone image LC. Besides, as was proved in Proposition [3l the charge-breaking 
stationary point is a minimum if and only if A'^'^ is positive definite in the entire Minkowski 
space, i.e. when all Aj < 0. 

If A^'^ is positive definite in the entire Minkowski space, the second matrix in fj4Tl) becomes 
positive definite, too. Therefore, in this case the neutral stationary point is a minimum if and 
only if > 0. But if C > 0, then lies outside the image of the lightcone LC. Indeed, 
rewriting ( IMl) with the lightcone vector = (1, ni, n2, n^) in the coordinate frame where 
Af^^ is diagonal, one gets 

^v\Ao, -\Ai\ni, -|A2|n2, -\A3\n3) - C(l, rii, ^2, ^3) = . 

In plain words, to arrive at starting from a point r'^ on LC, one must decrease its timelike 
coordinate Tq, while increasing the absolute values of all spacelike coordinates r'^. Such a shift 
brings us outside LC. Therefore, conditions for a charge-breaking minimum and neutral 
minimum are opposite and can never be met simultaneously. □ 

In other words, if one has a charge-breaking minimum, then the neutral stationary points 
are necessarily saddle points, and vice versa. This fact offers a very clear explanation to the 
finding that if a neutral minimum exists, the charge-breaking stationary points lie above it, 
see [13] and references therein. Note also that Proposition O implies that the charge-breaking 
or charge-restoring phase transitions cannot be of the first ordejf]. 

Another observation can be read off Eq. P2|) : 

Proposition 6. If the prototypical model has Mq < 0, then the neutral stationary point, when 
exists, is always a minimum. 

■^Here we mean the zero-temperature phase transitions from one minimum to another under the slow 
variation of the parameters of the Higgs potential. We do not consider here thermal phase transitions. 
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Proof. In deriving the mass matrix (1421) . we have aheady ahgned the timehke axis of A^^'^ with 
the future direction. To arrive at the prototypical model, one would need only a 3D rotation, 
which does not change Mq. Therefore, the value of Mq that enters (H2|) is the same as in the 
prototypical model. 

If Mo < 0, then ( = |f^Ao — Mq > 0, so that one needs to check only the eigenvalues 
of fH21) . But this matrix is a sum of a strictly positive and a non-negative matrices and, 
therefore, it is again positive definite. Its eigenvalues are all positive, and the stationary point 
is a minimum. □ 



4.2 Toy model 

The analysis of the number and properties of the solutions in the general case will be given 
later in geometric terms. Before we turn to it, it is instructive to consider in detail a toy 
model, which introduces these geometric constructions in a simple way. 

The toy model is a 2HDM potential with three equal eigenvalues of A^'^: Ai = A2 = A3 = A. 
Its prototypical model has A5 = and A = A3 + A4 and is invariant under any 3D rotation. We 
reiterate the point that the results obtained below hold not only for the prototypical model 
itself, but also for all the 2HDM's related to it by S0{1,3) transformations. 

The operator A^*^ in the toy model acting on the orbit space vectors preserves the 3D 
spherical symmetry. Consider the plane in the orbit space defined by the future 4-vector 
Cq = (1,0,0,0) and the mass terms 4-vectoil§ M^. The stationary points r^, if exist, must 
also lie on this plane. Therefore, following the strategy for the search of the stationary points 
outlined above, one must check only the two directions on the lightcone that belong also to 
this plane. 

The number of stationary points depends on the position of M'^, which is illustrated and 
explained in Figs. □ (A > 0) and [2] (A < 0). 

• If lies inside the future lightcone LC~^, then there are two neutral solutions. In terms 
of the components of M'^ in this particular basis, this occurs when rhfi, fh\2 > 0. 

• If lies outside LC~^ and outside the past lightcone LC~ (i.e. < implying 
and are of the opposite signs), then there is a unique solution. 

• If lies inside the past lightcone LC^ (^11 ? ^22 < 0); there are no solutions. Fig. [T]i. 

The saddle point/minimum classification of these points depends on the following three 
situations: 

• A > 0, which means A4 < and A3 = A — A4 > A; 

• A < and |A| < Aq, which means A > A3 > and A > A4 > 0; 

• A < but |A| > Aq, which means A4 > A and A3 < 0. 

Let us consider the first case, shown in Fig. [H If M^^ lies inside LC~^, we have two solutions, 
one on each of the two light-rays. Without lack of generality, we set mf^ > > (i.e. 

^If c>c Cq, there would be a whole sphere of solution to the stationary point problem. 
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Figure 1: Schematic representation of the search for the neutral stationary points for the toy 
model with A > 0. One takes a point r'^ on LC~^, finds its image r'^' on LC under action of 
A^'^ and then subtracts (n'^ to arrive at M^. Arrows going down correspond to C > 0, arrows 
going up correspond to C < 0. Shown are the four possible positions of M^: (a) M'^ inside 
LC, C > for both solutions; (b) inside LC but outside LC , C > only for solution 1; 
(c) M^" outside LC (M^ < 0), the only solution with C > 0; (d) inside LC", no solutions. 

M3 < 0). In the chosen basis, 

— 2 — 2 

TT) HI 

MX = Mo-M3 = ^, MX = Mo + M3 = ^, 

(n+|r2+) = Ao - A = A , = Aq + A = A3 , (e2|e2) = -A. 

The eigenvalues of the matrix (l42l) simplify to v'^{Aq — A) and twice degenerate value f ^(Aq — 
A) - 2Mo = 2C - v^A. The two solutions are 

2 M0-M3 (Mo-M3)2 
solution 1 V = 2— — , Vi = — 2- 



An-A ' Ao-A 



2 ^ Mo + M3 (Mo + M3)2 

solution 2 V = 2— — , V2 = —2 : . 43 

Ao-A ' Ao-A 

The first solution lies deeper than the second; in fact, the first solution is a minimum, since 

AMo - A0M3 . 

2C = 2 ^ ^ > and 2C - v^A = -2M3 > , 44 

Ao-A ' V / 

while the second solution is a saddle point due to 2( — v'^A = 2M3 < 0. 
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Figure 2: The same as in Fig. [T] but for A < 0, |A| < Aq. (a) M'^ inside LC", C < for both 
solutions; (b) inside LC but outside LC", C > only for solution 1; (c) outside LC 
(M^ < 0), the only solution with C > 0; (d) inside LC~, no solutions. 



If lies outside LC^, then fh\2 becomes negative, and the second solution disappears. 
The first solution, however, still corresponds to the minimum. 

Let us turn to the case A < 0, which means that A'"^ is positive definite in the entire 
Minkowski space. If |A| < Aq, then LC" still lies inside LC^ (this case is shown in Fig. [2]), 
while for |A| > Aq the cone LC" is wider than LC"*". 

The above expressions are still valid in this case, however the change of sign of A changes 
also geometric constructions. When lies both inside LC^ and LC", then, as before, two 
solutions exist, but both of them have C < 0, therefore, both correspond to the saddle points. 
In this case the charge-breaking stationary point is the only minimum of the potential, in 
agreement with Proposition [51 

For a neutral stationary point to become a minimum, C, must be positive, which happens 
when lies outside the lightcone image LC". Note that in the case the other neutral 
stationary point, if exists, must be a saddle point, since the mass matrices for these two points 
have among their eigenvalues 2M3 or — 2M3. 

It is interesting to trace how the properties of the vacuum change if A < is fixed and 
continuously shifts from inside to outside of the LC". Initially the ground state of the 
potential is charge-breaking, but upon crossing LC" it turns into a neutral one. The solution 
of the stationary point problem also moves continuously, and one observes a charge-restoring 
phase transition of the second order. We stress that this is a zero-temperature phase transition, 
not a thermal one. Such an evolution of M^, and therefore such a phase transition, might be 
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relevant for the early history of the Universe, |18j . 

We summarize the solution to the minimization problem in the toy model: 

• Non-trivial stationary points of the Higgs potential exist if and only if lies outside 
the past lightcone, in agreement with Proposition |H Among these stationary points, 
there is only one minimum. If other stationary points exist, they are saddle points. 

• If A > 0, the minimum of the Higgs potential is neutral. Depending on the position of 
M'^, it can be also accompanied by one neutral and one charge-breaking saddle points. 

• If A < 0, then the minimum can be either neutral or charge-breaking. If lies inside 
LC, the minimum is charge-breaking, while for outside LC, it is neutral. 



4.3 Geometric constructions in the toy model 

The above criteria for the existence and the number of solutions in the toy model can be 
obtained in a purely geometric fashion. 

As before, make appropriate boosts to set all Aqj = and consider the following con- 
struction. Take a sphere composed of points = if^(l,n) with fixed v"^ and all possible 
orientations of the unit 3D vector n. Under the action of A^'^ it is mapped to another sphere 
r''^ = |f^(Ao,An). At each point of this sphere attach now a straight line parallel to r^, i.e. 
consider the 3D manifold, defined by the l.h.s. of Eq.( 



Mv = Ix 



= (^Iv'Ao - C, (^Iv'A -C)^)] ^ (45) 

with v"^ fixed and all possible C ^ 1^ ^ind n. This manifold is a surface of a pair of oppositely 
oriented 4D cones, whose common vertex is located at ^(Aq — A, 0) lying on the future line. 

As f ^ grows, Aiy is homothetically rescaled. The shape of the cones remains the same, 
their common vertex sliding upwards. When f ^ changes from zero to infinity, each of the cones 
sweeps a region in the 4D space. The upper cones sweeps the interior of LC~^, while the lower 
cone sweeps all the space outside LC~ . Therefore, if a point is chosen inside LC^, it will be 
swept twice: once by each cone; a point outside LC~^ and LC~ will be swept only once, by 
the lower cone. This gives the number of solutions in the toy model. 

A slightly different way to arrive at this conclusion is to consider an "equal-time" 3D 
section of A/i^, that passes through the point M^, 



Mo-^t;2(Ao-A) 



n , (46) 



and see how it changes when f ^ grows from zero to infinity. If Mq is fixed and positive, is a 
sphere whose radius shrinks from Mq to zero and then swells all the way up to infinity. Points 
that lie inside the initial sphere (which correspond to lying inside LC^) will be swept 
twice; points outside that sphere will be swept only once. If Mq is negative, then the above 
section is a sphere whose radius starts from \Mq\ and increases to infinity. A point outside the 
initial sphere will be swept once, while a point inside it will never be crossed. 
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4.4 Geometric constructions in the general case 



We now turn to the general case with all Aj different. Without lack of generality, suppose that 
A3 > A2 > Ai. 

Again, consider the sphere = n) = |f^(l, Ui) for fixed v^. Upon action of A'^'^ it 

is mapped to the (surface of) ellipsoid |f^(Ao, Ajnj). The 3D manifold 



Mv = <x 



= Qt;^Ao - C, (^^'Ai - C) n^, (^^v'A, - n^, Qi^^As - c) ^^3) | (47) 

has a more involved geometry. It is composed of 4 parts, which are "glued" along 2D "equal 
time" regions: the upper cone-like manifold, which corresponds to ( E (— 00, |f^Ai], 

two compact cusped manifolds, A^l^^ and which correspond to C G [|f^Ai, |f^A2] and 

( G [|f ^A2, |f ^A3], respectively, and the lower cone-like manifold, Aii^^ corresponding to 

To help visualize this construction, we schematically show in Fig.[3]a similar warped surface 
in three-dimensional space, where there is only one compact manifold that links the two cones. 

The 2D regions along which these parts of Ai^ are glued, are ellipses together with their 
interiors. For example, the border between Ai^^ and M.^^ is an ellipse with semiaxes f ^(A2 — 
Ai) and f^(A3 — Ai) that lies in the plane orthogonal to the first axis. 




Figure 3: A schematic illustration of the lower-dimensional analogue of Ai^, ( H7I) . for some 
f ^ > 0. In this case it is a 2D ruled surface embedded into 3D space; it is made of the upper 
cone, the lower cone and one compact manifold that links them. The Mq = const sections are 
ellipses. 

Now consider the family of all Ai^ for = [0, 00). As f ^ changes, the above construction 
is rescaled. As before, the upper cone sweeps the interior of LC^ and the lower cone sweeps all 
the region outside LC~ . In addition, each of the two compact manifolds sweeps a certain cone 
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in the Minkowski space. These cones are the caustic 3D surfaces of the family of half-planes 
given by the l.h.s. of Eq. flMl) . In order to visualize these caustics cones and develop intuition, 
we suggest to take a look at a simple planimetric problem described in Appendix [Bl Returning 
then to 2HDM, consider again fj,^, the 3D section of M.y that passes through M^. It is an 
ellipsoid with semiaxes 

Mo-^v\Ao-A,), Mo-^v\Ao-A2), Mo - ^v\Ao - A^) . (48) 

If Mq is negative, the situation remains the same as before. As f ^ changes from zero to infinity, 
the ellipsoid defined by semiaxes fHHl) starts from sphere with radius \Mq\ and grows to infinity, 
sweeping exactly once each of the points that lie outside the initial sphere. Thus, if in 
the prototypical model lies inside the past lightcone LC~ , there are no non-trivial stationary 
points of the potential, in agreement with Proposition |H 

For Mo > there is a richer set of possibilities. Consider the three special values of v"^, 

l''i=J^.^ ^ = 1,2,3, 0<vl<vl<vl, 

when i-th semiaxis of fiy becomes zero. As long as < < w^, the 3D section passes through 
M^:\ It is an the ellipsoid that starts from the initial sphere, gradually shrinks, and at v"^ = vf 
it collapses to an ellipse with semiaxes 

_ A2 - Ai _ A3 - Ai 

that lies in a plane orthogonal to the first axis. Then, for vf < v'^ < the section passes 
through A/fl^^; the first semiaxis of grows while the other two shrink, and at = it 
collapses to another ellipse orthogonal to the second axis. Later, at = v^^ /i^, collapses to 
the third ellipse, and after that it keeps swelling up to infinity. 

The 3D section of the first caustic cone generated by A^i^'* is the region in swept by a 
family of ellipsoids with semiaxes (HHj) for f ^ < < w|. This region has a cusped astroid-like 
shapqj, somewhat elongated along the third axis, shown in Fig. HI The section of the other 
caustic cone has a similar shape, although oriented differently. 

One can show that as f ^ grows from vl to fg? each point inside the figure shown in Fig. H] 
is swept exactly twice. The same takes place for points inside the 3D section for the second 
caustic cone, when f ^ grows from f| to f |. In addition, each point in the 3D section considered 
will be always swept once for some v"^ > and, if ikf^ lies inside LC^, will be swept once 
more for some < v"^ < vf. 

The above construction gives a geometric proof of the following criteria: 

Proposition 7. (1) if M'^ lies outside LC~ , at least one neutral stationary point in the Higgs 
field space exists; 

(2) if lies inside LC^ , at least two neutral stationary points exist; 

(3) if M^, in addition, lies inside one of the caustic cones, two additional neutral stationary 
points appear; 

(4) if lies inside both caustic cones, four neutral stationary points appear, in addition to 
criteria (1) or (2). 



^In 2D case it would be simply a rescaled astroid, see Appendix [ 
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Figure 4: The Xq = const section of one of the cones of caustics. It is a 3D surface that 
envelops the family of ellipsoids fHSl) for vf < < v"^- 

Thus, the overall number of nontrivial stationary points can be as high as six neutral plus 
one charge-breaking, in agreement with results of ^12j . The sequence of these solutions, which 
are given by intersections of the evolving with a given point ikf^, is easy to visualize. For 
example, consider a generic that lies sufficiently close to the future-directed eigenvector 
of A'^'^. Then, it lies inside LC~^ and both caustic cones. As the ellipsoid starts from sphere 
and shrinks, it will pass this point once before collapsing to the first ellipse. Then, evolving 
between the first and the second mutually orthogonal ellipses, it will pass the point twice. 
Then, evolving from the second to the third ellipses it will yield two more intersections, and 
finally, swelling from the last ellipse to infinity it will pass the point the sixth time. 

A special discussion is needed when is not a generic point, but lies inside one of the 
ellipses described above, say, the first one. The planimetric problem described in Appendix [B] 
helps understand what happens in this case. Instead of two intersections for values of f ^ just 
below and just above vf, there are two different solutions at = vf, so that the counting 
rules stated in Proposition [7] remain unchanged. These two solutions correspond to non-zero 
and opposite values of the coordinate ri, breaking thus the Ci — > — Ci reflection symmetry of 
Af"" and M^. 

Note that such a breaking of discrete symmetry occurs if and only if the vector lies in- 
side at least one of the caustic cones. The spontaneous CP- violation in 2HDM to be discussed 
in the next section has precisely these geometric roots. 

Note that we analyzed above the most general case. Partially degenerate situations (e.g. 
Ai = A2, etc.) can be analyzed in a similar way. 

Proposition [7] gives the number of neutral stationary points of a given Higgs potential, but 
does not say anything about their minimum/saddle point nature. As was showed before, to 
establish that a given stationary point is a minimum, one needs to flnd all the eigenvalues 
of the mass matrix (I42p and make sure that all of them are positive. Whether this criterion 
admits a simple geometric formulation, remains to be seen. 
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5 Spontaneous CP- violation 



The Higgs potential is said to be explicitly CP conserving, if it commutes with the operator 
of the CP transformation. It is known that the Higgs potential in 2HDM is explicitly CP- 
conserving if and only if there exists such a SU{2) reparametrization transformation that 
makes all the coefficients real, see discussion in [8J. However, in this case it is possible that 
the vacuum expectation values for the two doublets attain a non-zero relative phase, which, in 
general, cannot be "rotated away" by a unitary transformation of the fields still keeping the 
potential real (for two examples of models when this phase can be rotated away, see fl9[ 120] ). 
When this happens, the CP-conservation is broken spontaneously. In fact, this possibility was 
the original motivation to consider EWSB with two Higgs doublets, [3]. 

Undoubtedly, for any explicitly CP-conserving Higgs potential with a given set of Aj and 
mjj there exists a unique answer to the question whether spontaneous CP violation occurs. 
However, up to now no such criterion for a general 2HDM has been published. What one 
encounters in literature is either straightforward minimization of the Higgs potential in some 
simple cases, [31 [191 EO], or discussion of what would be the properties of the CP-conserving 
or CP- violating vacuum solutions, if they exist, [H [11]. For example, [8] contains a set of 
basis-invariant conditions that must be satisfied in order to guarantee that no spontaneous 
CP- violation occurs. Aside from very complicated algebraic expressions for these conditions, 
they are presented in implicit form: in order to check if a given 2HDM exhibits spontaneous 
CP violation, one would need first to explicitly find the minima and then check if they satisfy 
the algebraic criteria. 

Here we tackle this problem with the geometric approach developed above. 

5.1 Preliminary conditions for spontaneous CP-violation 

Let us return for a moment to the conventional 5't/(2)-reparametrization group, which gen- 
erates 3D rotations in the orbit space. Here, instead of A^^, one deals with scalar Aqo, vector 
Aoi and rank- two symmetric tensor Ajj , [TOl [TT] . The four- vector is split into scalar Mq 
and 3D vector Mj. 

As shown in [TOj, the basis-invariant condition for the existence of a "hidden" real potential 
form (i.e. for the explicit CP-conservation), is that among the three eigenvectors of Aij there 
must be one orthogonal both to Aqj and Mj. We now extend this statement to the spontaneous 
CP- violation: 

Proposition 8. Spontaneous CP-violation can take place only if there exists one and only 
one eigenvector of Ay orthogonal both to Aqi and Mj . 

Proof. The first part is proved in \lOl. We now show that if there are two eigenvectors of Aij 
orthogonal both to Aqj and Mj, then spontaneous CP- violation cannot take place. 

The eigenvectors of a symmetric real matrix Aij are mutually orthogonal. If non-zero 
vectors Aoi and Mj happen to be orthogonal to two eigenvectors of Aij, then they must 
be proportional to the third eigenvector of Aij. In order words, Aqj and Mj are themselves 
eigenvectors of Aij. As shown in \iQ\, this situation corresponds precisely to the existence of 
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hidden Z2-symmetry in the Higgs potential. Diagonahzing Ajj, one arrives at 



/ An \ 

=\ A22 , Aoi = (0, 0, A03) , Mi = (0, 0, M3) . 

V A33 / 



We are now going to prove the following statement: even if one manages to find a solution 
with r2 = 2Im (0^02) 7^ 0, then one can always "rotate away" the relative phase still keeping 
the parameters real. Indeed, let us write the first and second components of the equation 



Two possibilities are to be considered: An = A22 and An 7^ A22- The first possibility implies 
also a possibility to have both rii and n2 simultaneously non-zero. However, the matrix Ajj 
possesses in this case an extra symmetry that allows one to manually align the first axis along 
the transverse vector n±. This rotation sets n2 = and removes the relative phase between 
the two v.e.v.'s. For the second possibility, one can have only one of ni and n2 nonzero, but 
not both. If this non-zero component happens to be n2, then rii is necessarily zero. One can 
then flip the first and second axes and arrive at the real Higgs potential with ^2 = 0. We 
conclude that there is no room for spontaneous CP-violation in this case. □ 

Put in simple terms, this Proposition means that in order to exhibit not explicit, but 
spontaneous CP- violation, the potential must have certain symmetry, but it must not be too 
symmetric. It might be interesting to see if this is a generic rule that is observed when a 
spontaneous breaking of a (discrete) symmetry takes place. 

The fact that the non-zero phase does not always imply spontaneous CP- violation has 
been known since long ago (see discussion in [5]). For example, authors of [TH] considered a 
2HDM with Ai = A2, Ae = A7 = real, A5 real, = 11122, ^12 ^^^l; showed that the 
relative phase between the v.e.v.'s can be rotated away still keeping all the parameters real. In 
our language, this model corresponds to = (Mi, 0, 0), Aoi = (Aqi, 0, 0) and diagonal Aj^, 
which obviously violates conditions of Proposition [SI Paper [20] deals with an even simpler 
example: = -^6 = -^7 = 0, which hterally corresponds to the situation just discussed in the 
proof of this Proposition. Our linear algebraic approach reveals the origin of these sporadic 
findings. 

5.2 Necessary and sufficient conditions for spontaneous CP- violation 

Given that the requirements stated in Proposition [H] are satisfied, the following questions now 
arise: (1) what are the necessary and sufficient conditions for existence of a spontaneous 
CP- violating stationary points and (2) when this stationary point is a minimum? 
The answer to the first question is given by the following Proposition: 

Proposition 9. Given that the conditions of Proposition are satisfied, the spontaneously 
CP-violating stationary points exist if and only if the components of the prototypical 
model satisfy the following inequality: 









Ml 



+ 



< 




(49) 
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Proof. As already shown at the end of the previous section, in order for the prototypical model 
to exhibit spontaneous CP-violation, must lie inside an appropriate caustic cone and on 
the plane orthogonal to the second axis. Such a non-trivial caustic cone exists if and only if 
A2 is a non-degenerate eigenvalue of A^'^. If this is satisfied, must lie inside the ellipse 
with semiaxes Mo|Ai - A2|/(Ao - A2) and M0IA3 - A2|/(Ao - A2), which gives (H9|l . 

To complete the proof, one needs to show that the passage from the original to the pro- 
totypical model and back respects the CP-symmetry properties. We start with the Higgs 
potential that has (hidden) explicit CP-conservation and bring it to the real potential form 
by an SU (2) reparametrization transformation. The corresponding 50(3) rotation in the orbit 
space brings the tensor A^" and 4-vector to the following form: 



/ • • 





■ \ 












-A2 





V ■ ■ 





■ / 



where dots represent generic values. Diagonalization of A^^" requires now only boosting along 
the first and third axes and 2D rotation between them. Such transformations do not generate 
the second components of Aq, or Mj, therefore the prototypical model is also explicitly CP- 
conserving. 

The inverse transformation laws, too, do not involve the second axis. If one obtains a 
(non)zero value of r2 in the prototypical model, this value remains the same upon returning 
to the original Higgs potential. Within the conditions of Proposition [HI the non-zero phase 
cannot be "rotated away", and the spontaneous CP- violation takes place in the original model 
if and only if it takes place in the prototypical model. □ 

In Appendix O we rederive condition ( H9l) with the aid of straightforward algebra. Unfor- 
tunately, it appears to be a non-trivial task to check if the conditions derived in [|8J for the 
spontaneous CP- violation in a generic explicitly CP-conserving model agree with ours. 

Appendix [C] lists also the necessary and sufficient conditions for the spontaneously CP- 
violating point to be a minimum. Here we cite only the particularly simple necessary condition: 
spontaneously CP-violating stationary point of the potential can be a minimum only if A2 is 
the largest spacelike eigenvalue of A^^ . 

When a CP- violating minimum is the global one still remains to be studied. 



6 Discussion and conclusions 

In this work we exploited the simple observation that the generic Higgs potential of the two- 
Higgs-doublet model keeps its form under any linear transformation between the two doublets. 
The most interesting subgroup SL{2, C) of this group of transformations induces the Lorentz 
group of transformations of = ^^cr^^. The Higgs potential takes a very compact and 
Lorentz invariant form V = —M^r^ + ^A^^^r'^r'^. This extended reparametrization symmetry 
turns out very useful in analysis of the existence and number of extrema of the Higgs potential 
and their classification according to neutral/charge-breaking and saddle/minimum nature. 

We introduced the notion of a prototypical model, which has only four free parameters in the 
quartic potential. All 2HDM's linked to the prototypical model by an SL{2, C) transformation. 
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share the same set of solutions to the stationary point problem. This construction gives the 
answer to the question which free parameters are physically essential and which are redundant 
when we study general properties of 2HDM. 

It has always seemed that explicit manipulation with high-order algebraic equations are 
unavoidable in order to learn something about the vacuum in 2HDM. The key point of the 
present work is the observation that we do not need such manipulations to understand the 
general structure of the 2HDM vacuum. 

Several geometric constructions appear naturally in the orbit space of 2HDM: the future 
lightcone LC'^ , its image under the action of A'^'^ , and the caustic cones. Using them, we 
proved that the charge-breaking and neutral minima cannot coexist in 2HDM. We also dis- 
covered the role of the caustic manifolds as separatrices between the 2HDM's with different 
number of solutions. Finally, we established the geometric roots of the spontaneous breaking 
of discrete symmetries and obtained explicit conditions when spontaneous CP- violation hap- 
pens in 2HDM. What remains to be seen is when these spontaneously violating minima are 
the global ones. 

A number of directions for further study opens up. 

• In the present work we did not consider in detail different types of degenerate situa- 
tions. It seems interesting to classify such situations in terms of little groups of residual 
symmetries, and trace down the consequences for the 2HDM vacuum. This study can 
be done within the geometric framework developed here, without any need of explicit 
expressions for positions of extrema. 

• If Aj are fixed and mfj continuously change, the global minimum undergoes certain zero- 
temperature phase transitions. Using geometric approach, one can classify the phase 
transitions allowed in 2HDM and study their properties. 

• Check if the geometric approach can be applied also to Higgs-Higgs scattering processes. 
In particular, search for a simpler rederivation of the perturbative unitarity constraints, 

• Study the deformation of the above description when quantum corrections to the poten- 
tial are taken into account. 

• Extend the geometric approach to other forms of the Higgs potential, for example, A^- 
Higgs doublet model. Indeed, in NHDM one can again construct a scalar and an N"^ — 1- 
component vector in the adjoint representation of SU (N), see pjj. Then, switching from 
reparametrization group SU{N) to SL{N,C), one can unite them in a single vector in 
1 + (A^^ — l)-space with Minkowski structure. This should provide extra aid to attempts 
to understand the general properties of NHDM, such as [13]. Finally, the formalism 
developed here might turn out useful for analysis of Landau-Ginzburg potentials in 
some condensed matter problems. 

I am thankful to Ilya Ginzburg, Roman Li, Celso Nishi, Alessandro Papa, Gianfranco 
Sartori, Gianpaolo Valente and the anonymous referee for discussions and useful comments. 
This work was supported by the INFN and FNRS, and partly by grants RFBR 05-02-16211 
and NSh-5362.2006.2. 
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A Diagonalization of 



We first show some basic facts on diagonalization of the real symmetric tensor A^'^ in the 
Minkowski space. A^*^ can be viewed as an operator acting on vectors in the Minkowski space. 
The (right) eigenvector with eigenvalue A is defined by 

A^Vp'^ = Ap^ 

thus the characteristic equation for the eigenvalues is det(A'^'^ — Kg^'^) = 0. One can switch to 
the more common expression by lowering one of the superscripts: 

det(A^^-A5^:) = 0. (50) 

Tensor A'^p that enters flSUl) is still real, but not symmetric anymore. Therefore, its eigenvalues 
are, in general, complex. If this is the case, the transformation that diagonalizes A^'^ would 
involve a complex transformation matrix, which does not belong to the S0{1,3) group. This 
means that exploiting only the proper Lorentz group of transformations, one might not be 
able to diagonalize a given A^'^. 

The situation simplifies due to the requirement that A'^'^ is positive definite on the future 
lightcon^: 

Proposition 10. Tensor A^^ is positive definite on the future lightcone if and only if the 
following three conditions hold: 

(1) A^^ is diagonalizahle by an 50(1,3) transformation, 

(2) the timelike eigenvalue Aq is positive, 

(3) all spacelike eigenvalues Aj are smaller than Aq. 

Proof. Obviously, if A^^ satisfies conditions (l)-(3), then the positive definiteness follows 
immediately. So, one needs to prove that conditions (l)-(3) do follow from the positive 
definiteness. 

Despite A^^ being real and symmetric, its eigenvalues can be complex because of the non- 
Euclidean metric. The first step is to prove that the positive definiteness in the future lightcone 
LC^ implies that all the eigenvalues of A^^j are real. 

Indeed, suppose there is a pair of complex eigenvalues, A and A*, with respective complex 
eigenvectors p^ and q^: 

A^^p^ = Ap^ , A^'^q" = \*q^ . 

One can show that there can be only one pair of complex eigenvalues, thus, A is non-degenerate. 
Since A(^ is real, q^ oc p'^* (and can be taken equal to p^*). These eigenvectors are orthogonal, 
p^q^ = 0, (it follows from the standard argument due to A 7^ A*) and can be normalized so 
that = g^g^ = 1. 

Consider now a real vector 

r'' = cp^ + c*p*^ . 

Suppose that r^r^ = + c*^ = 2|cpcos(20c) > 0, so that either or —r^ lies inside the 
forward lightcone. Then, the corresponding quadratic form is 

A^^r^r" = Ac^ + \*c*^ = 2|A| jcj^ cos (20^ + 0a) • 
am thankful to Celso Nishi who discovered a flaw in the earlier version of the proof of this proposition 
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Due to the phase shift 0a 7^ one can always find 0c such that cos(20c) > but cos(20c+0a) < 
0, i.e. one can always find an r'^ G LC^ such that K^^r^r^ < 0, which contradicts the 
assumption. 

Since all the eigenvalues of A^^, are real, the eigenvectors also can be chosen all real and 
orthogonal. One can show that they can be normalized so that one of the eigenvectors has 
positive norm, p^po^i = 1, while the other three have negative norms p^Pifi = — 1 for each 
2 = 1,2, 3. Thus, the transformation matrix T that diagonalizes K^y is real, and after di- 
agonalization A^j, takes form diag(Ao, — Ai, — A2, —A3). Note that transformation T also 
conserves norm, r^r^ = const. It means that T can be realized as a transformation from the 
proper Lorentz group. 

Now, the requirement that is positive definite in LC^ means 

Ao — p(Ai sin 9 cos + A2 sin 6^ sin + A3 cos > 

for allO<p<l,O<0<7r and all 0. This holds when Aq is positive and is larger than all 
Ai. □ 

Note that since A^*^ is a linear operator, it maps points inside some surface S into points 
that also lie inside the image of this surface 5*'. 

The positivity conditions in the form of inequalities among A's can be readily obtained 
from this statement. Here we reproduce the standard set of positivity conditions for a model 
with explicit Z2-symmetry, i.e. when Ag = A7 



0. of this model has form 



/ ^ + A3 



2 



v - 










A4 + Re A5 — Im A5 
— Im A5 A4 — Re A5 

- A3 y 



which after a rotation among the first and second axes and an appropriate boost along the 
third axis gives 



/ A3 + yiA^ 



\ 





A4 + As 

A4 - A5 

- A3 / 









Such a transformation is possible if and only if Ai + A2 > |Ai — A2I, which is equivalent to 
Ai > 0, A2 > 0. The eigenvalues of A^^, are 



2An 



A3 + V A1A2 , 2Ai = — A4 — IA5 



2Ao 



-A4 + A5 



2A3 



A3 — ^/ A1A2 



The conditions stated in Proposition [TO] now translate into 

Ai>0, A2>0, A3 + \/AiA2 > 0, A3 + A4 - IA5I + VA1A2 > 0, 
which coincides with the standard set of inequalities. 
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B Planimetric primer: caustics, number of solutions 
and spontaneous violation of discrete symmetries 

The following simple planimetric problem provides insight into phenomena that take place in 
the original problem. 

Given a linear map of the 2D vector space rj — > r ■ = AijVj — Bi with real symmetric matrix 
Aij and real -Bj, how many points are there on the unit circle such that r' is proportional to 

t K 1 




Figure 5: (Left) The caustic curve formed by the family of straight lines AijTj — C-ri for lying 
on the unit circle. A point in the region bounded by this curve has four lines passing through 
it. For a point Bi deliberately chosen on the horizontal axis, there are nontrivial solutions, 
one of which is also shown, that do not lie on it. (Right) The same solution visualized as a 
shrinking ellipse. When an ellipse collapse onto the line, on which B^ lies, two points from 
AijTj arrive at Bi. 

Following the strategy outlined in Section 14.11 take a point on the unit circle and draw 
a straight line AijTj — ( ■ Vi trying to arrive at Bi. Repeating this for every point on the unit 
circle, one gets a family of straight lines. If Aij is not proportional to the unit matrix, these 
lines do not intersect at a single point, but instead form a caustic curve (an astroid in the 
present case). If Bi is inside the region bounded by the caustic curve, there are four solutions; 
if Bi is outside the caustics curve, there are only two solutions. Fig. [5l left, illustrates this 
construction. For a point inside the astroid deliberately chosen to lie on the horizontal axis, 
there are two "trivial" solutions, which also lie on this axis, and two nontrivial solutions, one 
of which is also shown. 

Although both Aij and Bi are symmetric under flip of the vertical axis, the two non-trivial 
solutions are not. Thus, having Bi inside the caustic curve makes spontaneous breaking of a 
discrete symmetry possible. 

In order to understand better the origin of this phenomenon, consider another way to 
visualize the solution of this problem. The matrix Aij acting on the circle rj maps it to an 
ellipse with semiaxes ai and a2, the eigenvalues of Aij. As ( increases, the ellipse AijVj — ( ■ Vi 
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shrinks, see Fig. 0, right. At C = Oi and C = ct2 this elhpse collapses to an interval on the 
vertical and horizontal axes, respectively. In the latter case, in particular, one obtains -82 = 
without requiring r2 to be zero. 



C Spontaneous CP- violation in the prototypical model: 
straightforward algebra 



Here we obtain the criterion ( l49i) using straightforward algebra similar to the original calcu- 
lations of [3]. Suppose A5 > A4, which guarantees that the non-trivial minima of the potential 
are neutral and allows us to search the solutions in form 

^.-^=J1.A- (51) 



V2\ pi J ' V2\ P2e 

After substituting them into the Higgs potential of the prototypical model (fT5l) . we obtain 
(suppressing barred notation for simplicity) 

V = -^{ml^pl + ml^pl + 2ml^piP2Cos9) + ^ [X{p^^ + pt) + 2 plpj^X^ + X4 + X5 cos 29)] . (52) 

This potential can have both CP-conserving and CP-violating stationary points. Let us 
find conditions when spontaneous CP-violation takes place. Since sin 6' 7^ 0, we first set 
dV/dcosO = 0, which gives 

2 



The condition for non-zero phase is therefore 



\6\<1 - ^<4p?p^ (54) 



2 \2 
12_ 
2 
5 



The potential can be then rewritten as 



1 1 1 (m^ 

V = ~{ml,pl + ml,pl) + [X{p\ + p\) + 2pIpIX,,,] + ^A5P?p^(cos ^ - 5f - , (55) 

where A345 = A3 + A4 — A5. Upon differentiating the potential in respect to p\ and P2 (note 
that both pi and p2 are non-zero) and by using cos 9 — 5 = 0, one obtains: 



2 _ "m/V — '"-22^345 2 _ "^22^ ~ "^11^345 /p-^x 

Pi - 72 ' P2 - 72 ■ \^^) 



2\ 2\ 2\ 2\ 

2 _ ^11^ ~ '"^22^345 ^2 _ ""^22^ ~ "^11^345 

P 

Exploiting the identity 



X2 \2 ' r-z ^2 \2 

^ ~ ^345 ~ ^345 



4:{ax — by){ay — hx) = (a — 6)^(x + yY — (a + 6)^(x — yY , 
one rewrites (15^ as 

{ml2? {mil + "^^2? _ (^11 - "^22)^ /.^x 
Ai (A + A345)2 (A-A345)2 ■ ^ ^ 
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This coincides with 
Ao 



if one recalls the definition of and notes that 



A2 = ^(A + A345), A2-A3 = i(A-A345), Ag - Ai = A5 • 



In order to establish the minimum conditions for the CP-violating stationary point, we cal- 
culate the matrix of second derivatives of the potential in respect to pi, p\ and 9 



( 



D 



A 

A345 



A345 
A 



— sin6'm?a^ \ 
12 pi 



sin6'm^2~ 



P2 



(5J 



— sin6'm^2^ ~sin6'mf2^ 4 sin^ 6'p^p2A5 y 
This matrix is positive definite if and only if 

A > , A5 > , A^ - A^45 > , det D > . 



(59) 



From the first three conditions one obtains the following necessary conditions for the sponta- 
neous CP- violating point to be a minimum: A2 > Ai and A2 > A3. The last condition in fISUl) 
can be rewritten as 

("^12) ^ ^ „2 „2 A 



< 4ptP2 



^345 



A, 



A(4 + 4) -2A 



345 



(60) 



compare it with flM|) . This extra condition cuts out a subregion in the ellipse defined by 
where the spontaneous CP- violating solutions are minima. Using the parameters of the 
prototypical model, one can define 



X 



Ml A0-A2 
Mo ' A2 - Ai 



M3 A0-A2 



Mo Aa-A, 



The condition (JSlj) for the CP- violating stationary points is then + < 1, while the 
condition fl60p for the CP- violating minimum takes the form 



A2-A1 ^ A2-A1 
Ac - A2 A2 - A3 



< (1- z 



2\2 



(61) 



Fig. [6] illustrates the shape of this region in a representative case Aq = 4, Ai = 1, A2 = 3, 
A3 = 0. 

Note that checking whether the CP- violating minimum is a global one is still a non-trivial 
task even in the prototypical model. 
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